Gravitation model for spatial network based on the heterogeneous node 
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In this paper we consider nodes in network are heterogeneous and the link between nodes is caused 
by the potential dynamical demand of the nodes. Such demand can be measured by gravitation 
which increases with the heterogeneous strength of node and decreases with the geographical dis- 
tance. Based on this, we propose a new model for spatial network from the view of gravitation. The 
model is to maximize the potential dynamical demand of the whole network, indicating the possible 
maximal efficiency of the network and the highest profits that operators may gain. The model can 
vary its topology by changing two parameters. A simulation for the Chinese city airline network is 
completed. In the end of this article we discuss the significance and advantage of the heterogeneous 
nodes. 

PACS numbers: 89.75.Hc, 89.75.Da, 89.40.Dd 



I. INTRODUCTION 



II. HETEROGENEOUS NODES AND 
GRAVITATION 
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Since the initial studies on the small-world phe- 
nomenon were presented by Watts and Strogatz [l| and 
the scale- free property by Barabasi and Albert , a lot of 
achievements on complex network have been gotten. And 
our research group have also studied some features on 
network Most previous works focus on the topolog- 
ical properties of the network. However many networks 
are those embedded in the real space whose nodes occupy 
a precise position in Euclidean space and whose links are 
constrained by the geographic distance. Some typical ex- 
amples_are communication networks 0, Q , electric power 
grids transportation systems ranging from river Q 
to airport d, H, street [ll|, railway and subway 
In these spatial networks, geographical factor is 
demonstrated to play an important role on the network's 
topology. Very recently a few models for the spatial 
network have been proposed. Some models consider 
both the topological and spatially preferential attach- 
ment 1131 1 1 511 while others take some optimal mech- 
anism [la , Il7l . [l8| . These works provide some guidelines 
in network design. On the other hand, the geographic 
effects on the efficiency or traffic of the network are paid 
less attention. However traffic on network is very im- 
portant since it represents the efficiency of the network. 
Therefore, if the traffic between nodes can be predicted in 
some way, it is likely to construct the network efficiently. 
Inspired by this idea and its significance, we proposed a 
new spatial network model. The model is to maximize 
the whole expected traffic of the network, indicating the 
highest efficiency that the network may gain. The ex- 
pected traffic is measured by the gravitation. Especially, 
we propose the heterogeneity of nodes and argue its sig- 
nificance on the network's topology and dynamics. 
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In the previous studies on complex networks, nodes are 
usually identical. However in the real network, situations 
are very different. In the city airline network, each node 
(city) has different population or economic level and such 
difference can be very obvious. As is shown in Fig. [TJ 
the population of the Chinese city follows approximately 
a two-regime power-law distribution, indicating a large 
variance of the grade of the city. The population can be 
considered as the own attribute of node which is distin- 
guished from others and indicates its own grade. Such 
phenomenon also exists in the Internet and World Wide 
Web-routers with different capacity and webs with het- 
erogeneous resources. Thus a node usually has its own 
attribute, grade or energy indicating the heterogeneity. 
We define heterogeneous strength to describe this het- 
erogeneity, denoted by M. In addition, we found the het- 
erogeneous strength(i.e.the population of city node) have 
positive correlation with the degree of city node, as is 
shown in Fig. [2] 

Now we consider a network is composed of a group of 
nodes and links. Links connect every node and realize 
the dynamics among them. If two nodes with no dy- 
namics between them, the link of the two vertices can be 
thought unwanted. In a circuit network, for example, a 
lead equals to disconnection if its current is 0. In other 
words, it is the dynamical demand between nodes that 
causes the link. Such demand is fulfilled when nodes are 
connected and can be measured by the weight. However, 
the demand still exists even though the vertices are not 
connected. 

Consider a network with n nodes and edges. Every 
two nodes have their demand for some information ex- 
change. What we care is which of these demands are 
the most exigent. If such information is obtained in 
advance, the network will be constructed efficiently by 
preferentially investing those pairs of nodes with great 
demand. To measure such demand, some of its proper- 
ties should be discussed first: (1) In a spatial network, 
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can be described by the following equation, 
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FIG. 1: The Zipf-type plot distribution of the population of Chi- 
nese cities. It approximately follows a two-regime power-law distri- 
bution indicating the large variance of the heterogeneous strength. 
The population information is obtained from Ref. [23l 



Gi 



MfMf 
K — - — — 

D 1 - 



(1) 



Equation (1) reminds us of the Newton's gravitation 
equation. From the gravitation view, the dynamical de- 
mands is considered to make the nodes magnetize with 
each other and get the link when the gravitation Gij is 
great enough. However such demand is potential, it is a 
measurement of the expected traffic which is a prediction 
of the real traffic. 



III. GRAVITATION MODEL FOR SPATIAL 
NETWORK 

We use the gravitation equation (1) to measure the 
potential dynamical demand. Rewrite equation (1), 
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FIG. 2: The correlation of the degree and the population of the 
city in Chinese city airline network. It shows the heterogeneous 
strength has a positive correlation with the node's degree. The de- 
gree information is obtained from URL: [http: //w ww. caac.gov.cn 



links are constrained by the geographic factors and the 
travel cost among the nodes increases with the distance 
[HI]. So for the spatial network, the demand of nodes 
is assumed to decrease with the distance between them. 
(2)Since such demand can be measured by weight after 
nodes are connected, we hope it has some similar features 
with the real weight. In some kinds of network, such as 
the airline network, the weight has the form as follows, 
< Wij >^ Xij(kikj) a (l9l . [20 |. where fej and kj are the 
degree of node i and j. Since degree k has a positive cor- 
relation to heterogeneous strength M as we mentioned 
above, it is reasonable to suppose such demand is pro- 
portional to (MiMj) a . The above two points indicate 
that the potential dynamical demand between two nodes 
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Mj are the het- 



where K is a constant coefficient. Mi 
erogeneous strength of node i and node j. Exponents a, 
7 determine the impact of Dij and M on the gravita- 
tion. If 7 takes great value, the network is expected to 
be strongly constrained by the geographic distance. On 
the other hand, larger value of a means the dominant 
impact of the heterogeneous strength on the topology of 
the network. The topology of the simulated network is 
changed by varying parameters a and 7. The cost of 
constructing the network is defined as the total number 
of the network's edges. It is reasonable for the airline 
networks whose cost and expense are related to the to- 
tal hops. But for road network whose expense depends 
on the Euclidean distance, such definition seems uncon- 
vinced. However, later discussion will be revealed that 
such definition is still reasonable. Now suppose there 
are n nodes distributed on a two-dimension plane. The 
heterogeneous strength and coordinates of each node are 
known. By equation (2), the gravitation Gij of any two 
nodes can be calculated. We connected preferentially 
those pairs of nodes with greater gravitation Gij and 
completed such process when the link comes to the value 
we preset (namely the cost). Such process can be de- 
scribed as an optimal model, 
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where fcy is adjacency matrix element of the network. 
The upper part of the equation means to take the maxi- 
mum and the bottom part means the maximum process 
with that constraint. However, the above process may 



FIG. 3: Four networks with different topology are shown which 
is controlled by the parameters a and 7. (a) a takes any value 
and while 7=0. The network is dominated by two large hubs when 
7=0; (b) a = 1 and 7=1. With the increasing of 7, the hubs 
become much smaller than figure (a); (c) a = 1 and 7 = 2; (d) a 
= 0, 7 takes any value. The topologies are strongly reminiscent of 
airlines and roads, respectively. 

cause some isolated nodes. So two more restrictions are 
introduced to ensure each node is connected, 

Y^hj > 10' = 1,2, 3, ..,n) 

i 

Y, k H >l(* = l,2,3,...,n) (4) 

i 

Following the method above, a network of thirty nodes 
is simulated. The thirty nodes are distributed on the two- 
dimension plane randomly and each node is assigned a 
heterogeneous strength. Set the cost e = 39 and the 
coefficient K = I (actually K makes no difference to the 
model). By varying the value of a and/or 7, we got four 
networks with different topology as is seen in Fig. [3] 

For Fig. 3(a), the value of 7 = makes the network 
only rely on the heterogeneous strength of nodes, This 
causes the topology is dominated by two hubs since the 
nodes of great heterogeneous strength is easy to magne- 
tize others. Whereas with the 7 increasing, more effect 
of the geographic factor makes the node tend to con- 
nect to the closer ones and weaken the hub-and -spoke 
effect. When a—1 and 7=2 (Fig. 3(c)), the network ex- 
hibits some features similar to the airline network. When 
a—0 (Fig. 3(d)), the topology is entirely constrained by 
the geography, which forms a two-dimensional network 
strongly reminiscent of roads. As we have defined in the 
section II, the gravitation Gy describes some dynamical 
demands or expected traffic. From the view of opera- 
tor, the significance of the equation (3) is to maximize 
the whole expected traffic of the network which reflects 
the efficiency of the network. In other words, our model 
provides a possible way for operators to gain the highest 
profit. Even though a link costs much, the operators will 



FIG. 4: The geographic distribution of the cities with airports. 
Each city denotes a node. There are totally 121 nodes in our sim- 
ulation. 

still earn the cost and gain profit in short order as long 
as the traffic within the link is large enough. Thus such 
a link will be constructed all the time. In summary, the 
cost is a minor factor compared to the dynamical demand 
and that is the reason we simplify the cost. 



IV. SIMULATION FOR THE CHINESE CITY 
AIRLINE NETWORK 

The gravitation model we proposed is used to simu- 
late the Chinese airline network [2l| . The cities are con- 
sidered as the nodes of the network while the edges re- 
flect the airlines. To reproduce the real network better, 
we set the number of nodes n=121 and number of links 
E=689. The position of each node is sited on the two- 
dimensional plane, as it is shown in Fig. U] The heteroge- 
neous strength M and the distance Dij are, respectively, 
defined as the population of the city and the Euclidean 
distance of city i and j. Selecting a—1 and 7 €[1,2] [22| . 
they are to describe the interactions of cities. Here we 
set 7 = 1.5. 

Fig. [5] shows the simulated network. Obviously, the 
hubs in the real network such as Beijing, Shanghai, 
Guangzhou, Harbin and Urumchi exhibit the similar hub- 
and-spoke phenomenon in our simulation. In real con- 
dition, Beijing, Shanghai and Guangzhou are the three 
cities with the highest degree while in our model they are, 
respectively, Shanghai, Beijing and Wuhan (Guangzhou 
is the fourth). The reason for such difference may be 
that using the population to denote the heterogeneous 
strength of city is intuitive but not be exact because the 
economy and the administration factors are also impor- 
tant indexes for the grade of city. In spite of this dif- 
ference, we still succeed in reproducing the every hub 
and their hub-and-spoken phenomenon existing in the 
real network. The average shortest-path length L and 
degree-degree correlation exponent r of the model net- 
work are calculated, where L = 2.302, r = -0.401 while 
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FIG. 5: The simulation for the Chinese airline network. Each link 
of two nodes reflects an airline between them. There is totally 689 
links. The cities, such as Beijing, Shanghai, Guangzhou, Harbin 
and Urumchi, exhibit hub-and-spoke phenomenon in our simulation 
just as they do in the real network. 
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FIG. 6: The clustering-degree distribution of the simulation. The 
distribution satisfies a linear decreasing feature with slope -1 in 
log-log coordinate, which indicates that the result of simulation 
reproduces the hierarchy. 



FIG. 7: The degree distribution of the simulation. It presents 
a behavior of two-regime power-law degree distribution with the 
exponent7i = -0.46 for the first power laws and 72 = -2.3 for the 
second. The turning point happens at degree k = 22. In compari- 
son, 71 = -0.53, 72 = -2.05 and the turning point at degree k = 20 
in the Chinese City Airline network. Both the exponents and the 
turning point fit the real network well. 



strength of our simulation. Considering the gravitation 
is an effective prediction of the real traffic, we let the 
gravitation of the two nodes represent the weight of the 
link between them and the node strength is the sum of 
such gravitation among all the links to it. As is shown in 
Fig. [51 the node strength increases with the degree, but 
does quicker than linearly, as 1.45 power, namely satis- 
fies S(k) ~ k 1A5 . This result reflected the real network 
well. In the Austrian airline network (24l|and in the Chi- 
nese city airline network it indicates that node strength 
increases with the degree in a non-linearly manner. 



MORE DISCUSSIONS ABOUT THE 
HETEROGENEOUS NODES 



in the real network L = 2.263, r = -0.408. Fig. [5] shows 
the clustering-degree distribution. It meets C(k) ~ fc" 1 
which indicates the model network exhibits the same hi- 
erarchy as the real network does. 

Fig.[7]is the degree distribution of the model network.lt 
satisfies the two-regime power-law distribution. This is 
an interesting result not only to verify the real degree dis- 
tribution but also to produce the power-law distribution 
without any preferential attachment mechanism. We are 
not sure if it is a coincidence or caused by some other in- 
ternal reasons, but we will still try to offer an explanation 
in the following discussion. 

Fig. [5] shows the correlation of the degree and node 



As we have discussed in section II, the nodes of many 
real networks are usually heterogeneous. But one may 
ask if it is necessary or meaningful to take this factor into 
consideration. We argue the heterogeneity of nodes may 
be the key reason causing the topology that we observed 
in many networks. Some clues are gotten from Fig.[TjThe 
population of the city obeys a two-regime power-law dis- 
tribution which is similar to the degree distribution of 
the city airline network. This reminds us if it is the 
heterogeneous nodes that cause the heterogeneous net- 
works whose degree distribution has a large variance or 
in other words, homogeneous nodes will hardly produce 
the hubs or scale-free properties. An indirect evidence is 
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FIG. 8: The correlation of the degree and node strength follows 
S(k) ~ fc 1,45 , indicating the non-linearly correlation. 



work's evolvement. Now it may be understood why our 
model can produce the power-law distribution without 
the preferential attachment mechanism. Because such 
mechanism has been produced by gravitation although it 
is not introduced designedly. Constrained by the geog- 
raphy, nodes with small heterogeneous strength usually 
connect to the large ones nearby while the latter are able 
to link with other nodes of even larger heterogeneous 
strength far away and so on. It can be imagined that 
the network will exhibit the hierarchy by such organi- 
zational manner. Besides heterogeneity of nodes is also 
related to the dynamics - the traffic. As equation (2) 
shows, the heterogeneous strength influences the gravi- 
tation dj which is considered as the expected traffic of 
node i and j . Although it is only a prediction of the real 
traffic, Fig. [5] shows that such prediction verifies the real 
condition. This point indicates that the heterogeneous 
strength of the node is possible a bridge by which we can 
associate the topology with the dynamics. 

VI. CONCLUSION 



that hubs in the real network are usually those with large 
heterogeneous strength. If admitting the heterogeneity of 
nodes, some phenomena such as hubs and preferential at- 
tachment can be explained from the view of gravitation. 
Considering a node with large heterogeneous strength in 
a spatial network, it produces a larger gravitation field ac- 
cording to the gravitation view. For its extensive area of 
influence, many nodes are magnetized and connected to 
it. So the nodes with large heterogeneous strength usu- 
ally have larger degree, which indicates the emergence of 
hubs. When a new vertex emerges, it is more possible 
to be magnetized by the nodes with large heterogeneous 
strength in the same way, which makes it tend to con- 
nect to these nodes (as is mentioned above, these nodes 
usually have larger degree). And this just causes the phe- 
nomenon of the preferential attachment observed in net- 



In this paper, link is caused by the potential dynamical 
demand among nodes which makes the nodes magnetize 
with each other and its value can be estimated by gravita- 
tion. Based on this, a new model for the spatial network 
is proposed. It can generate the topology of airline and 
road network and can behave well in simulating the real 
Chinese airline network. The gravitation model has its 
practical meaning that follows a principle of maximiz- 
ing some expected benefit to construct and optimize the 
network. We argue that cost is a minor factor despite 
its impact of restriction compared with the dynamical 
demand of nodes. Heterogeneous nodes and gravitation 
view are two ideas indicated by our model. We think 
they are important in understanding the topology and 
the dynamics of the network. However, more demonstra- 
tion studies are essential to support this point of view. 
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